Geometrical exponents of contour loops on synthetic multifractal rough surfaces: 

multiplicative hierarchical cascade p model 
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In this paper, we study many geometrical properties of contour loops to characterize the mor- 
phology of synthetic multifractal rough surfaces, which are generated by multiplicative hierarchical 
cascading processes. To this end, two different classes of multifractal rough surfaces are numerically 
simulated. As the rst group, singular measure multifractal rough surfaces are generated by using 
the p model. The smoothened multifractal rough surface then is simulated by convolving the rst 
group with a so-called Hurst exponent, H* . The generalized multifractal dimension of isoheight 
lines (contours), D{q), correlation exponent of contours, xi , cumulative distributions of areas, ^, 
and perimeters, rj, axe calculated for both synthetic multifractal rough surfaces. Our results show 
that for both mentioned classes, hyperscaling relations for contour loops are the same as that of 
monofractal systems. In contrast to singular measure multifractal rough surfaces, H* plays a lead- 
ing role in smoothened multifractal rough surfaces. All computed geometrical exponents for the 
rst class depend not only on its Hurst exponent but also on the set of p values. But in spite of 
multifractal nature of smoothened surfaces (second class), the corresponding geometrical exponents 
are controlled by H* , the same as what happens for monofractal rough surfaces. 



I. INTRODUCTION 

Random phenomena in nature generate ubiquitously 
fractal structures which show self-similar or self-aSine 
properties [l|-y| . When the fractal structure of a system 
is uniform and free of irregularities, we have monofractal 
structure. A monofractal system can be characterized 
by a single scaling law with one scaling exponent in all 
scales. For a self-affine surface and interface, this expo- 
nent is called roughness exponent or Hurst exponent, H . 
Surface with larger H seems locally smoother than the 
surface with smaller if @, S] ■ 

In topics rangin g fr om biology [5|, |6|, surface sciences 
THlOl . turbulence pH - flSj . diffusion-limited aggregation 
14| . bacterial colony growth [i5|, climate indicators [ia| 
to cosmology [17[ , there are many surfaces and interfaces 
exhibit multifractal structures. A multifractal system 
can be considered as a combination of many different 
monofractal subsets [3, Q . Multifractality manifests it- 
self in systems with different scaling properties in various 
regions of the system. In addition, multifractals can be 
described by infinite different numbers of scaling expo- 
nents h(q), where q can be a real number. The appear- 
ance of the infinite different numbers ensures that the 
theoretical and the numerical study of the multifractal 
surfaces is more complicated than those of monofractal 
ones. Changing one of the h{q)'s can lead to different 
feature in the system. One of the important character- 
istics of the multifractality is the presence of the sin- 
gularity spectrum, /(a), which associates the Husdorff 
dimension /(a) to the subset of the support of the mea- 
sure fj, where the Holder exponent is a; in other words 



/(a) = diniH {x\fJ.{Bx{e)) ^ e'*}, where Bx{e) is an e-box 
centered at x. 

A single scaling exponent can be determined for a 
monofractal structure, by use of various methods [1, [a, 
[l8l - [2a | . Not only a spectrum of exponents but also differ- 
ent algorithms should be computed for a multifractal fea- 
ture (power spectral, distribution method and so on) and 
these may give different results for a typical multifractal 
case [2J]. Thus, the better and more complete theoreti- 
cal frame work, the better our understanding, providing 
deeper insight to observational multifractal rough sur- 
faces. 

Recently, isoheight contour lines has been utilized to 
explore the topography of rough surfaces and it exhib- 



ited interesting capabilities 2J-|31|. The contour plot 
consists of closed non-intersecting lines in the plane that 
connects points of equal heights. The fractal properties 
of the contour loops of the rough surfaces can be de- 
scribed by just the Hurst exponent [30, [3l|- This result 
was confirmed in different systems with quite different 
structures in recent years both experimentally and nu- 
merically. Using numerical approach, the predicted rela- 
tions were confirmed in glassy interfaces and turbulence 
[32| , in two-Dimensional fractional Brownian nrotion |26| , 
in KPZ surfaces |25| and in discrete scale-invariant rough 
surfaces [27[. The predictions were also confirmed by us- 
ing experimental data coming from the AFM analysis of 
WO (3) surfaces [H. 

However, although there have been many studies con- 
cerning the contour lines of monofractal rough surfaces, 
there ares neither theoretical nor numerical inferences 
about the contour lines of multifractal rough surfaces. 



Because of the presence of numerous exponents in the 
multifractal surfaces, theoretical study of multifractal 
surfaces seems to be difficult. Moreover, in many previ- 
ous methods, the exponents determined by fractal anal- 
ysis, generally provide information about the average 
global properties, whereas geometrical analysis addresses 
information from point to point. J. Kondev et al. 
pointed out that geometrical characteristics can discrim- 
inate various monofractal rough surfaces that have a sim- 
ilar power spectrum [3l|, [S^ • Therefore, the geometrical 
properties may introduce a new opportunity to charac- 
terize multifractal surface. 

It is worth noting that because contour sets are the 
intersection of a horizontal surface in particular height 
fluctuation and do not reflect the full properties of fluc- 
tuations in various scales, it is not trivial that the geo- 
metrical properties of multifractal rough surfaces based 
on isoheight nonintersecting feature behave in a multi- 
fractal manner as well. Therefore we use a new approach 
to investigate these processes. 

In this paper, we try to investigate multifractal struc- 
tures utilizing contour loops. We study the multifractal 
properties of a particular kind of multifractal surfaces. 
Two different type of synthetic multifractal rough sur- 
faces, namely singular measure and smoothened features, 
are generated. Two mentioned types have a multifrac- 
tal nature. Despite the complexity nature of the model, 
hyperscaling relation is satisfied for both categories. In 
addition, from contouring analysis, all of geometrical ex- 
ponents for various smoothened multifractal rough sur- 
faces are controlled by corresponding so-called Hurst ex- 
ponent, H*, This is also is the same as what happens for 
mono-fractal cases. However, for a singular measure mul- 
tifractal rough surface, geometrical exponents depend on 
the set of pvalues that is used to generate the underlying 
rough surface based on multiplicative cascade model. 

The structure of the paper is as follows: In the next 
section we will review multifractal rough surfaces. The 
hierarchical model to generate the surfaces will also given 
in this section. The multifractal detrended fluctuation 
analysis in two dimensions that is used to characterize the 
multifractal properties of rough surfaces will be explained 
in section HI. In section IV, nonlinear scaling exponents 
of multifractal rough surfaces are introduced. Section 
V will be devoted to numerical results for determining 
the scaling exponents of the contour loops in multifractal 
rough surfaces. In the last section we will summarize our 
findings. 



II. MULTIFRACTAL ROUGH SURFACE 
SYNTHESIS 

Recently there has been an increasing interest in the 
notion of multifractality because of its extensive appli- 
cations in different areas such as complex systems, in- 
dustrial and natural phenomena. Dozens of methods 
for the synthesis of multifractal measures or multifrac- 
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FIG. 1: Upper panel: Different steps of generating multifrac- 
tal rough surface in one Dimension. Lower panel: The same 
steps for multifractal rough surface in two diinensions [11]. 



tal rough surfaces have been invented. One of the most 
common methods thath can be followed deterministically 
and stochastically is the multiplicative cascading pro- 
cess [1, [nl, [22I, IsJ- Some of these synthesis methods 
are known as the random /? model [i2[, ol model j35j . 
log-stable models, log-infinitely divisible cascade models 
[361 137| and p model [11| . They were successfully applied 
in the studies related to rain in one dimension, clouds 
in two dimensions and landscapes in three dimensions as 
well as many other fields J36l - l39l | . 

The p model method was proposed to mimic the kinetic 



energy dissipation field in fully developed turbulence [11 



The so-called p model represents the spatial version of 
weighted curdling feature and is known as conservative 



cascade . It is based on Richardson's picture of energy 
transfer from cores to fine scales base on splitting eddies 
in a random way [40| In this model there is no divergency 
in corresponding moments in contrast to the so-called 
hyperbolic of a model [111, |4l[ 

On the other hand, many scaling exponents of men- 
tioned model can be determined analytically; therefore, 
it is a proper method to simulate synthetic multifractal 
processes ranging from surface sciences and astronomy 
to high energy physics such as cosmology and particle 
physics e.g QCD parton shower cascades, and cosmic mi- 
crowave background radiation J43 - |44| . In the context of 
p model simulating a synthetic one-dimensional data set, 
consider an interval with size L. Divide L into two parts 
with equal lengths. The value of the left half corresponds 
to the fraction < p < 1 of a typical measure \x while 
the right hand segment is associated to the remaining 
fraction (1 — p). By increasing the resolution to 2~", the 
multiplicative process divides the population in each part 
in the same way (see the upper panel of Fig. [T|). 

To simulate a mock multifractal rough surface in two 
dimensions, one can follow the same procedure as above. 
Starting from a square, one breaks it into four sub- 
squares of the same sizes. The associated measures for 
each cell at this step are pi/z for the upper right cell, •p2[i. 
for the upper left cell, ps/x for the lower right cell and p4/i 
for the lower left cell. The conservation of probability at 
each cascade step is pi +P2 +P3 +P4 = 1- This partition- 
ing and redistribution process repeat and we obtain after 
many generations, say n, 2" x 2" cells of size IjL = 2~" 
(see lower panel of Fig. [J). In the stochastic approach, 
the fraction of measure for each sub-cell at an arbitrary 
generation is determined by a random variable A with a 
definite probability distribution function P{A). By redis- 
tribution of measure, based on independent realization of 
the random A at smaller scales, one can generate a ran- 
dom singular measure over a substrate with size L x L 
as 



1(1) 



Hn{r;l) =Ain A(r) 



n{l) = log2 - ^ oo, (1) 



where r shows the coordinate of the underlying cell with 
size I. In this work, we rely on the stochastic version 
of the cascade p model to generate the synthetic two- 
dimensional multifractal rough surface (see Figs. [2] and 
[3]). The probability distribution function for our ap- 
proach is given by 



The so-called multifractal scaling exponent, T{q), and the 
generalized Hurst exponent, h{q), are quantities that rep- 
resent the multifractal behaviors of rough surfaces (see 
section III for more details). For p model cascade, these 
exponents can be calculated explicitly. The scaling ex- 
ponent T{q) is defined via partition function as 



n(l) 



ZS) 



limV|P(A,Or~r('?^. 



;-!-o 



(4) 



Using the value of P{A), e.g. for binomial cascade model 
P{A) = \[5{A~p) + 5{A-{l-p)], one finds 



T{q) = 



lim 

(^0 



log(^.(0) 
log(0 



{E-l){q-l)-\og^{p'i + {l-pY), (5) 



where E is the dimension of the geometric support, where 
for our rough surfaces is E' = 2. For generalized p model, 
the analytic expression of multifractal scaling exponent 
in two dimensions is given by [45| 



r{q) = -\og2{p\+Pl+Pl+Pl)- 



(6) 



One can use the above theoretical expression to get the 
most relevant quantities of the multifractal behavior and 
check the reliability and the robustness of numerical 
method. 

Recently, factorial moments, G-moments, correlation 
integrals, void probabilities, combinants and wavelet cor- 
relations have been used to examine maiiy interesting fea- 
ture of multiplicative cascade processes [4a] . But there is 
some ambiguity in properties of such processes that rep- 
resent multifractal phenomena. On the other hand, sen- 
sitivity and accuracy of results are method dependent; 
consequently, it is highly proposed to simultaneously use 
various tools in order to ensure the reliability of given 
results for underlying multifractal rough features. More- 
over, to make a relation between experimental data and 
simulation, generally, we require more than one charac- 
terization [ail m 

In the next section, to investigate the multifractal 
properties of simulated rough surfaces in two dimensions, 
we will introduce the so-called multifractal detrended 
fluctuation analysis. 



III. MULTIFRACTALITY OF SYNTHESIS 
ROUGH SURFACE 



P{A) = -[S{A-Ai) + S{A-A2) 
+6{A-A3) + S{A-A4)], 



where 



^1 

^3 
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P3, 



A2 

Aa 
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There are many different methods to determine the 
multiscaling properties of real as well as synthetic niul- 

(2) tifractal surfaces such as spectral analysis ^g, fluctua- 
tion analysis |49| . detrended fluctuation analysis (DFA) 

ISOl Isii. wavelet transform module maxima (WTMM) 
Il3,ll3, m m and discrete wavelets [H, IHj. For 
real data sets and in the presence of noise, the multifrac- 

(3) tal DFA (MF-DFA) algorithm gives very reliable results 




FIG. 2: Left: Contour plot at some typical levels of a singular 
multifractal rough surface generated by binomial cascade mul- 
tifractal method with p = 0.22(77 = 0.803). The right panel 
indicates the contour lines of the same surface convolved with 
H* = 0.700. The system size is 256 x 256. 



pM ISOJ I . Since it does not require the modulus maxima 
procedure therefore this method is simpler than WTMM, 
however, it involves a bit more effort in programming. 

In this work, we rely on the two dimensional multi- 
fractal dctrended fluctuation analysis (MF-DFA) to de- 
termine the spectrum of the generalized Hurst exponent, 
h{q). We then compare given results with theoretical 
prediction to check the reliability of our simulation. Sup- 
pose that for a rough surface in two dimensions height 
of the fluctuations is represented by T-L{y) at coordinate 
r = (i, j) with resolution A. The MF-DFA in two dimen- 
sions has the following steps [34| 

Stepl : Consider a two dimensional array 'H{i,j) where 
i = 1,2,...,M and j = 1,2,..., TV. Divide the H(i,j) into 
Ms X Ns non-overlapping square segments of equal sizes 
s X s, where M^ = [-^] and Ns = [^]. Each square 
segment can be denoted by Hi^^w such that Hi/,^(i,j) = 
'H{li + i,l2 + j) for 1 < i,j < s, where li ^ {v — l)s and 
h = [w - l)s. 

Step 2: For each non-overlapping segment, the cumu- 
lative sum is calculated by: 



Yu^wihj) = X! X! ^!^/!i'(^i'^2); 



(7) 



fei=l k2 = l 



where 1 < i,j < s. 

Step 3: Calculating the local trend for each segments 
by a least-squares of the profile, linear, quadratic or 
higher order polynomials can be used in the fitting pro- 
cedure as follows: 



Bu,wii,j) = ai + bj + c, 



(8) 
(9) 



Then determine the variance for each segment as follows: 
T^u,w{i,i) = Y^^^{i,j) ~ By^.u,ii,j), (10) 

^'»(^) = ^EE^'-(*'^')- (11) 

A comparison of the results for different orders of DFA 
allows one to estimate the type of the polynomial trends 
in the surface data. 





FIG. 3: Upper panel: A part of height fluctuations of singular 
measure mentioned in Fig. 2. Lower panel: The same surface 
convolved with H* = 0.700. 



Step 4 ■ Averaging over all segments to obtain the gr'th 
order fluctuation function 



F,is) = 



1 



Ms Ns 

Ms X TV, ^ ^^ 



;^'»(s)] 



9/2 



1/9 



(12) 



where Fq(s) depends on scale s for different values of q. 
It is easy to sec that Fq{s) increases with increasing s. 
Notice that Fg{s) depends on the order q. In principle, q 
can take any real value except zero. For q ~ Eq. p^ 
becomes 



Fo(s) = exp 



1 



Ms Ns 



2M, X N, 



j:j:^nFUs) 



(13) 



For (7 = 2 the standard DFA in two dimensions will be 
retrieved. 



Step 5: Finally, investigate the scaling behavior of the 
fluctuation functions by analyzing log- log plots of Fg{s) 
versus s for each value of q, 

i^(s) -s''(9'. (14) 

The Hurst exponent is given by 

H = h{q = 2)-1. (15) 

Using standard niultifractal formalism [50| we have 



Tig) = qhiq) - E. 



(16) 



It has been shown that for very large scales, N/A < s, 
Fq{s) becomes statistically unreliable because the num- 
ber of segments Ng for the averaging procedure in step 4 
becomes very small [3J]. Thus, scales N/A < s should be 



excluded from the fitting procedure of determining h{q). 
On the other hand one should be careful also about sys- 
tematic deviations from the scaling behavior in Eq. (fT2|) 
that can occur for the small scales s < 10. 

The singularity spectrum, /(a), of a multifractal rough 
surface is given by the Legendre transformation of T(q) 
as 



/(a) = qa-T{q), 



(17) 



where a = a ■ It is well-known that for a multifractal 

dq 

surface, various parts of the feature are characterized by 
different values of a, causing a set of Holder exponents 
instead of a single a. The interval of Holder spectrum, 
a S [amin,Q!max], Can be determined by J56ll57l| 



a„ 



lim 



= lim 



dT{q) 



driq) 



g->--oo dq 



(18) 
(19) 



To evaluate the statistical errors due to numerical cal- 
culations we introduce posterior probability distribution 
function in terms of likelihood analysis. To this end, 
suppose the measurements and model parameters to be 
assigned by {X} and {6}, respectively.The conditional 
probability of the model parameters for a given observa- 
tion is as follows (posterior) 



p{e\x) 



c{x\e)p{e) 
J c{x\e)p{&)d&' 



(20) 



here C{X\Q) and P(8) are called Likelihood and prior 
distribution, respectively. The prior distribution contain- 
ing all initial constraints regarding model parameters. 
Based on the central limit theorem, Likelihood function 
can be given by a product of gaussian functions as fol- 
lows: 



ln£(X|e) 



-x'(e) 



(21) 




FIG. 4; (Color online) Diagrams of h{q) (upper panel) and 
T{q) (middle panel) for different surfaces. We have distin- 
guished different surfaces with their Hi = hi{q = 2) — l coming 
from the Table I. The subindex (i £ [1, 12]) of each Hi (Hurst 
exponent) throughout this paper corresponds to a given set 
of p values reported in Table I. The lower panel corresponds 
to the singularity spectrum of a typical multifractal rough 
surface with H4 — 0.608. In all diagrams, symbols and solid 
lines correspond to results given by numerical calculation and 
theoretical formula, respectively. 



where e.g., for determining h{q) we have {X} : {Fq{s)} 



as observations and {0} 
determined. Also 



{h{q)} as free parameter to be 



X'{h{q))=j2 



[i^obs.(s)-Fthe.(s;/i(g))]^ 



^obs. 



is) 



(22) 



where Fobs.(s) is computed by Eqs. p^ and (|13p . 
FthcX^j h{q)) is the fluctuation functions given by Eq. 
([T4|. The observational error is crobs.(s)- By using the 
fisher matrix, one can evaluate the value of the error-bar 
at la confidence interval of h{q) [5^ 
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Ha 


= 0.706 
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0.100 


Hu) 


= 0.697 


0.120 


0.180 


0.560 


0.140 


i/11 


= 0.806 


0.160 


0.180 


0.170 


0.490 


Hu 


= 0.906 


0.410 


0.200 


0.210 


0.180 



TABLE I: The p values used for construction of surfaces with 
various Hurst exponents, Hi = hi{q = 2) — 1. The subindex 
(i € [1, 12]) of Hi represents the label of different sets of p 
values. 
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'Tiq) 



VW) 



(24) 



Finally we report the best value of the scaling expo- 
nent at la confidence interval according to h{q) ± a{q). 
Using the method mentioned in the previous section, we 
simulated multifractal rough surfaces and checked their 
multifractality nature by using the spectrum of h{q). Fig- 
ure|4]shows the generalized Hurst exponent and T{q) as a 
function of q for various values of measure sets reported 
in Table H The subindex (i £ [1, 12]) of each Hi (Hurst 
exponent) throughout this paper corresponds to a given 
set of p values reported in Table HI In addition the singu- 
larity spectrum of a typical simulated multifractal rough 
surface has been shown in the lower panel of Fig. ID 
The q dependence of h{q) as well as the extended range 
of singularity spectrum demonstrate the multifractality 
nature of synthesis rough surfaces. Theoretical predic- 
tions of t((7), /i(g) and f{a) shown by the solid lines in 
the corresponding plots, are given by Eqs. (HI), (|T6)) and 
(fT7|) . respectively. There is a good consistency between 
theoretical predictions and computational values. 

Before going further it is worth mentioning that in the 
cascade p model for various sets of p values which have 
the same h{q = 2), in principle, there exist different h{q) 
spectrums. To show this point, we fixed the value of 
T{q = 2) in Eq. ^ and by having e.g., pi and p2j one 
can compute the rest of p values according to normaliza- 
tion of p's. In Fig. [5] we show MF-DFA results of various 
sets of p values causing the same so-called /i(g = 2) ex- 
ponents. Subsequently, it is expected that for character- 
izing the geometrical properties of underlying surfaces, 
one must take into account full spectrum of generalized 
Hurst exponents. 

It must be pointed out that the generated surfaces have 









p =0.090, p,=0.100, P3=0.610, p,=0.200, H,=0.608 
p =0.600, p;=0.100, P3=0.237, Pj=0.063, Hj=0.608 
Pi=0.350, p;=0.100, P3=0.546, p^=0.004, H,=0.608 
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Pj=0.350, Pj=0.100, P3=0.546, Pj=0.004, H,=0. 
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FIG. 5: (Color online) The multifractal spectrum of surfaces 
produced by different sets of p values but with the same h[q = 
2) up to our numerical precision. 



some discontinuities (see Fig. [2]). To make them smooth, 
a proper way is using fractionally integrated singular cas- 
cade (FISC) method ^M,- In this method, the multifrac- 
tal measure is transformed into a smoother multifractal 
rough surface by filtering the singular multifractal mea- 
sure [/i(r), (Eq. ([T]))] in the Fourier space as 



H{v) = ^(r) ® |r| 



-{l-H') 



(25) 



where ® is the convolution operator and H* E (0, 1) is 
the order of smoothness (sec the right panel of Fig. [5] 
and lower panel of Fig. [3]). In this case Tf{q) reads as 



vii) = •^(9) + qH* 



(26) 



where T(q) is given by Eq.®. Using the correlation func- 
tion, C(|r|) ^ |r|~''', and its Fourier transform one can 
derive the power spectrum scaling exponent /3 of the sin- 
gular as well as the smoothened synthetic multifractal 
surfaces. To this end we demand the scaling behavior for 
power spectrum to be 



S{k) 



where k 



run from 1 to A^ = 



i^y } ■! i^x 



(27) 
2^^ j and (i,j) 



AxAf"' ^y ~ AxN 

L/A (the pixel of system size). Sub- 
sequently the power spectrum scaling exponent is given 
byiH 



P 



l + 2F*-l0g2(p2 + (l_p)2). 



(28) 
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FIG. 6: (Color online) Upper panel: Pprofile of singular (left) 
and smoothened (right) multifractal rough surfaces along a 
typical horizontal cut in Fig. 2. Lower panel: Spectral density 
of mentioned mock rough surfaces. The solid lines in the lower 
panel corresponds to a power law fitting function and symbols 
are given by numerical calculation. Here we took H* = 0.700. 



To make more sense, in Table HIl we collected the corre- 
lation and power spectrum exponents of stochastic pro- 
cesses in one and two dimensions . 



Exponent 
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ID-fBm 
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2D-fBm 
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2-2H 


-2H 


1-2H 


-1-2H 
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2H - 1 


2H+1 


2H 


2H + 2 



TABLE II: The most relevant exponents concerning stochas- 
tic processes in one and two dimensions. 

Figure |6] indicates one dimension profiles obtained 
along a typical horizontal cut in Fig. [2] for singular and 
smoothened multifractal rough surfaces. The lower pan- 
els of Fig. inishow the power spectrums of simulated rough 
surfaces. The convolution does not change the multifrac- 
tality nature of singular measure (see Fig. [7|). In this 
plot one can see that the synthetic smoothened surface 
remains multifractal. 



IV. GEOMETRICAL EXPONENTS OF 
CONTOUR LOOPS 

For a given multifractal rough surface with the height 
'H(x), a level set 7^(x) = "Ho for different values of Hq 
consists of many closed non-intersecting loops. These 




FIG. 7: (Color online) Generalized Hurst exponent of sin- 
gular measure for Hq = 0.802 (square symbols) and that of 
convolved with H* = 0.700 (circle symbols). The solid lines 
are from the theory. 



loops are recognized as contour loops. The contour loop 
ensemble corresponds to contour loops of various level 
sets. In Fig. [5] we plotted a set of contour loop at some 
typical levels for singular multifractal rough surface and 
corresponding convolved surface with H* = 0.700. The 
loop length s can be defined as the total number of unit 
cells constructing a contour loop multiplied by lattice 
constant A. The radius of a typical loop is represented 
by R and it is the side of the smallest box that com- 
pletely enwraps the loop. For a mono-fractal surface, 
these loops are usually fractal and their size distribution 
is characterized by a few scaling functions and scaling 
exponents. For example the contour line properties can 
be described by the loop correlation function G{r). The 
loop correlation function measures the probability that 
the two points separated by the distance r in the plane 
lie on the same contour. Rotational invariance of the 
contour lines forces G{r) to depend only on r = |r|. This 
function for the contours on the lattice with grid size A 
and in the limit r 3> A has the scaling behavior 



G(r) 



-2X1 



(29) 



where xi is the l oop correlation exponent. It was shown 
numerically [271 ISOl . |31| that for all the known mono- 
fractal rough surfaces this exponent is superuniversal and 



equal to i. A key consequence of this result is that, 
the contour loops with perimeter s and radius R of such 
surfaces are self-similar. When these lines are scale in- 
variant, one can determine the fractal dimension as the 
exponent in the perimeter-radius relation. The relation 
between contour length s and its radius of gyration R is 



{s){R)^R' 



(30) 



where Df is the fractal dimension and R is defined by 



i?2 



N Z-^i=l \^' 



=)' + {v. 



Vc 



with 



-N 1^1=1 ^i '"^"^^ Vc = W J2i=i Vi being the central mass 
coordinates. The Df is the fractal dimension of one 
contour and for mono fractal rough surfaces is given by 
Df = ^^^ [3l| . Depending on what feature of the mul- 
tifractal rough surface is under investigation one can get 
various types of fractal dimensions. In this paper we in- 
troduce the fractal dimension of a isoheight line, Df , and 
the fractal dimension of all the level set, d. The general- 
ized form of fractal dimension can be expressed by means 
of partition function of underlaying feature, which is con- 
tours in this context, as 



^^' i^aq-l \og{l) ■ 



(31) 



where I is the size of the cells that one uses to cover the 
domain and its minimum value is equal to grid size, A. 
Zq{V) is the partition function defined in Eq. (j?]) but here 
it should be constructed by using contour loops instead of 
height function and q can be any real number. It is easy 
to show that Di^q = 0) = D f and Di^q = 1) corresponds 
to the so-called entropy of underlying system [45| . 

For a given self-similar loop ensemble, one can define 
the probability distribution of contour lengths P{s). This 
function is a measure for the total loops with length s and 
follows the power law 



P{s) 



(32) 



where 77 is a scaling exponent. Another interesting quan- 
tity with the scaling property is the cumulative distribu- 
tion of the number of contours with area greater than A 
which has the following form 



P>(A)~A- 



(33) 



For mono fractal rough surfaces we have ^ = 2 — H. Using 
the scaling property of the mono-fractal surfaces it was 
shown that the three exponents Df, r?, ^ and xi satisfy 
the following hyperscaling relations [30[ 



D 



2x1 -2 



77-3 



(34) 
(35) 



Using the above relations it is easy to get the relation 
between ry and Hurst exponent H. Before closing this 



Exponent 


Relation 


Description 


^i 


G{r) ~ r-2-1 


Loop correlation exponent 


Df 


{s){R)^R^f 


Fractal dimension of a contour loop 


D{q) 


Eq. (31) 


Multifractal dimension 


d 


N{1) r^l"' 


Fractal dimension of all contour set 


V 


P{s) ~ s-" 


Length distribution exponent 


e 


P^{A) ~^-3 


Area cumulative exponent 



TABLE III: The relevant exponents introduced in this paper 
to characterize synthetic multifractal rough surfaces. 



section, we summarize all of the exponents introduced in 
this section in Table Hill 

In the next section we will calculate all mentioned ex- 
ponents by using different numerical methods for singular 
as well as smoothened multifractal rough surfaces and we 
will examine the validity of the hyperscaling relations in 
this context. 



V. NUMERICAL RESULTS 

In order to examine the geometrical exponents of the 
contour loops mentioned in Table IIIII of synthetic mul- 
tifractal rough surfaces, we have generated multifractal 
rough surfaces with different h(q = 2)'s using the typ- 
ical measures reported in Table HI We have generated 
100 ensembles of each surfaces with various sizes rang- 
ing from (2048 x 2048) to (4096 x 4096). To extract the 
contour loops of the inock multifractal rough surfaces at 
mean height. Ho, we use two different methods, the con- 
touring algorithm and Hoshen-Kopelman algorithm [2al . 
According to our results, these two methods give almost 
the same results for geometrical exponents. In the next 
subsections we present our numerical results concerning 
the exponents introduced in the preceding sections. 



A. Loop correlation function exponent 

The loop correlation function exponent xi is the most 
central exponent in mono-fractal rough surfaces. It is 
independent of H and is equal to ^ . This result has also 
been proven for iJ = according to the exact solvable 
statistical mechanics model for contours equivalent to the 
critical 0(2) loop model on the honeycomb lattice [3l|, 
[53. 

To find the correlation function from a given loop en- 
semble for multifractal rough surfaces, we followed the 
algorithm described in Ref. [3i|. We calculated the loop 
correlation function G{r) for our multifractal rough sur- 
faces (with system size 2048 x 2048 and averaging is done 
over 10 realizations). The log-log diagram of G{r)r'^^' 
versus r for different sets of p values {Hi = hi{q = 2) — 1 
for some i G [1, 12]) have been shown in the Fig. [5] Each 
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FIG. 8: (Color online) Log-log diagram of r'^^'G{r) versus r 
for different Hurst exponents. Upper panel corresponds to 
singular measure with the sets of p values reported in Ta- 
ble I. Lower panel is indicates loop correlation function for 
smoothened multifractal surface for various _ff*'s. In these 
figures we shifted the y axis vertically. The system size is 
4096 X 4096. 



set corresponds to a synthetic multifractal rough surface 
generated according to the algorithm presented in sec- 
tion IL Our results demonstrate that the xi exponent, 
not only depends on the value of the Hurst exponent but 
also depend on the different sets of p values (see the up- 
per panel of Fig. |8] ). In other words, as reported in 
Table U as well as shown in upper panel of Fig. (|S]), the 
sets i = 4, i = 5 and z = 6 of p values have equal Hurst 
exponent, nevertheless, the corresponding correlation ex- 



FIG. 9: (Color online) The log-log plot of {s){R) versus R for 
synthetic multifractal singular rough surface for H4, = 0.608. 



ponents, xi, for these sets differ completely. On the other 
hand, at the level of our numerical accuracy, as shown in 
the lower panel of Fig. [51 the value for the smoothened 
multifractal surfaces correlation exponents is the same as 
that reported for the mono fractal rough surfaces, namely 
xi = i 



B. Fractal dimension 

To calculate the fractal dimension of a contour loop, 
we have calculated the perimeter and radius of gyra- 
tions of different contour loops. Figure [9] shows log-log 
plot of {s){R) versus (R) values for synthetic multifrac- 
tal rough surfaces with typical value of Hurst exponent, 
H4 = 0.608. There are two distinct regions with differ- 
ent slopes in the diagram; the first region is related to 
a large number of small loops with radius smaller than 
one {R < 1) with £>/ = 1.00 ± 0.01. This is not a rele- 
vant phenomenon and it comes from the contouring al- 
gorithm that produces lots of contour loops around very 
small clusters (made usually from one cell). In the sec- 
ond region {R > 1) the slope increases to 1.43 ±0.02 and 
it keeps to follow the scaling behavior up to very large 
sizes. The slopes for different Hurst exponent follow the 
relation Dj = (3 — H)/2 for mono-fractal case [31|. For 
various values of the Hurst exponent our computation is 
shown in the upper panel of Fig. [TU] (see also Table ITV]) . 
At 1(7 confidence interval all slopes are the same. On the 
contrary, in the case of the contour lines of the convolved 
rough surfaces with arbitrary H*'s the fractal dimension 
of a contour line follows the formula of a mono fractal 
surface with H ^ H* , namely Dy = (3 - H*)/2. It is 
quite interesting that these results are completely inde- 
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FIG. 10: (Color online) Upper panel: Log- log of {s){R) versus 
R for singular multifractal rough surfaces for various sets of p 
values reported in Table I. Lower panel: The same diagram for 
smoothened synthetic multifractal rough surfaces. The sam- 
ple size is 4096 x 4096 and the ensemble average was done over 
100 realizations. To make more sense, we shifted the values 
of (s) vertically for different multifractal rough surfaces. 



pendent form the p values (lower panel of Fig. llOp . This 
simply means that the fractal dimension of the contour 
loops of the singular rough surfaces does not change with 
respect to the h{q = 2). In other words, in contrast to the 
mono fractal case, h{q = 2) alone can not represent the 
properties of the underlying singular multifractal rough 
surface. 

Wc also calculated the fractal dimension by using par- 
tition function introduced in Eqs. (|1]) and (PT|) . Fig- 



FIG. 11: (Color online) Generalized fractal dimension versus 
q for singular measure with H4 = 0.608 and that of convolved 
by H* = 0.700. For singular and smoothened surfaces Df = 
1.46 ± 0.05 and Df = 1.19 ± 0.05, respectively. 



ure [TT] shows D{q) as a function of q. The q depen- 
dance of these results confirms that contour loops of syn- 
thetic singular and smooth multifractal rough surfaces 
are multifractal. For q = at 68% confidence interval 
D{q = 0) = 1.46 ± 0.05. This is also in agreement with 
the value determined by calculating the scaling behavior 
of the contour sizes. In addition as wc may expect this 
diagram demonstrates that the isoheight contour loops of 
underlying simulated multifractal rough surfaces behave 
as a multifractal feature . 

As mentioned , the fractal dimension of all the con- 
tours, d, differs from the fractal dimension of a contour 
loop Df. The fractal dimension of a contour set for 
mono fractal rough surfaces is given hy d = 2 — H. For 
the smoothened multifractal rough surfaces introduced 
by Eq. (1^51). t he fractal dimension of the contour set is 
d = 2 — H* [13[. We have calculated the fractal dimension 
of the contour set by using the box counting method. As 
previously, we used a least-squares equation (Eq. (|21[) ) to 
determine the slope in the log-log diagram of the number 
of segments that will cover the underlying feature N{1) 
versus length scale I for different Hurst exponents. To 
obtain best fit value for the slope corresponding to our 
data, as well as its error, we divided the data into dif- 
ferent ranges and determined the slope by least-squares 
method. To do so according to likelihood function (Eq. 



T])), we define x as 



M 



x'{d) = Y. 



[N{h)~Nti,,.{k;d)f 



i=l 



aihy 



(36) 



where Af is the number of partitioning, namely J\f 
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FIG. 12: (Color online) Fractal dimension of all contours of 
the smoothened multifractal surfaces as a function of H* . The 
solid line corresponds to linear fitting function. 



L/l_\f, Nthc.ili',d) ^ l~ and a^{li) is the variance of the 
data in the corresponding range. Finally we determined 
the niinimum x^ ^-nd the best slope for the data. Fig- 
ure [12] corresponds to synthetic smoothened multifractal 
rough surfaces. In addition we checked that whether the 
result associated to the smoothened rough surfaces de- 
pends on the set of p values correspond to the same value 
of h{q = 2). Our findings confirm that d doesn't depend 
on different sets of p values. However, for the singular 
measure, d depends on the value of H and even p's used 
for the cascade algorithm. It has no regular behavior with 
respect to h{q = 2). Moreover for various sets of p values 
giving the same value of h{q ~ 2), one finds out different 
values for fractal dimension of all contour sets. This is 
quite surprising because for singular measure multifrac- 
tal surface, we have H* = and, therefore, if the formula 
2 — H* was correct in this regime, we should have d = 2 
for all the different h{q — 2)'s. We are not aware of any 
theoretical argument that can explain this phenomenon. 
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FIG. 13: (Color online) Upper panel: The cumulative distri- 
butions of the areas of the contour loops with respect to the 
area for the singular multifractal rough surfaces. The corre- 
sponding set of p values is given in Table I. Lower panel: The 
same distribution for the smoothened multifractal surfaces. 
For clarity, we shifted the value of y axis vertically for both 
diagrams. 



To calculate the exponent ^ we have calculated the 
P>(A)A^/^ with respect to the area of the contour loops. 
In Fig. [T3]and table II VI we have shown the results for vari- 
ous values of Hurst exponent reported in table H] and aver- 
aging is done over 100 realizations. The results are quite 
different from what we expect for mono- fractal rough sur- 
faces. For mono fractal rough surfaces we have ^ = 2~H. 
It must be pointed out that for synthetic singular multi- 
fractal rough surfaces ^ decreases by increasing H, which 
is the same as mono-fractal rough surfaces. In addition ^ 



not only depends on h{q = 2) but also is affected by other 
values of /i(g)'s. This finding is due to the multifractality 
nature of the singular measure rough surface. The same 
computation for the smoothened multifractal rough sur- 
faces is shown in lower panel of Fig. [131 This results 
confirm that the exponent is controlled by H*, and ^ is 
given by the same equation as for the mono fractal rough 
surfaces. 
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FIG. 14: (Color online) Upper panebThe perimeter distri- 
bution exponent for different sets of p values of the singular 
measure. Lower panel the same measure for the smoothened 
multifractal rough surfaces. The values of y axis are shifted 
vertically. 



D. Probability distribution of contour length 



H 


V 


Df 


C 


2x1 


Hi = 0.305 


2.67 ±0.03 


1.43 ±0.04 


2.44 ± 0.06 


1.60 ±0.10 


H2 = 0.404 


2.60 ±0.03 


1.41 ±0.04 


2.30 ± 0.06 


1.49 ±0.10 


H3 = 0.504 


2.50 ±0.03 


1.42 ± 0.04 


2.16 ±0.06 


1.25 ±0.05 


H4 = 0.608 


2.45 ±0.02 


1.42 ± 0.04 


2.04 ±0.06 


1.30 ±0.03 


H5 = 0.608 


2.74 ±0.02 


1.43 ±0.04 


2.50 ±0.06 


1.70 ±0.03 


He = 0.608 


2.64 ±0.02 


1.42 ± 0.04 


2.31 ±0.06 


1.53 ±0.03 


Hs = 0.706 


2.35 ±0.02 


1.43 ± 0.04 


1.90 ± 0.06 


1.12 ±0.03 


Ha = 0.802 


2.27 ±0.02 


1.44 ±0.04 


1.80 ± 0.06 


1.02 ±0.03 



TABLE IV: Different geometrical exponents of the contour 
loops extracted from surfaces with different sets of p— values 
reported in Table I. Theses values completely dependent on 
the p values. 



H 


7,-1 




3Df + 2x1 


DfV + 2 


Hi = 0.305 


1.67 ±0.03 


1.71 ±0.06 


5.89 ±0.16 


5.82 ±0.12 


H2 = 0.404 


1.60 ±0.03 


1.63 ±0.06 


5.72 ±0.16 


5.67 ±0.12 


H3 = 0.504 


1.50 ±0.03 


1.52 ±0.06 


5.51 ±0.13 


5.55 ±0.11 


Hi = 0.608 


1.45 ±0.02 


1.44 ±0.06 


5.56 ±0.12 


5.48 ±0.10 


H5 = 0.608 


1.74 ±0.02 


1.75 ± 0.06 


5.99 ±0.12 


5.92 ±0.11 


Hg = 0.608 


1.64 ±0.02 


1.63 ± 0.06 


5.79 ±0.12 


5.75 ±0.11 


Hs = 0.706 


1.35 ±0.02 


1.33 ± 0.06 


5.41 ±0.12 


5.36 ±0.10 


Hg = 0.802 


1.27 ±0.02 


1.25 ±0.05 


5.34 ±0.12 


5.27 ±0.10 



TABLE V: Verification of two basic hyperscaling relations for 
synthetic singular measure multifractal rough surfaces. 



cal exponents of the contour loops of mono-fractal rough 
surfaces and singular multifractal rough surfaces, the hy- 



perscahng relations 



Df 



77- 



I and D 



f 



2X1 



V-3 



- are valid 



up to numerical accuracy (see Table. IIVI and Table. \V\ . 
The important factors in obtaining this hyperscaling re- 
lation concern power- law relations for P{s) and P>(A). 
The second hyperscaling relation comes from the follow- 
ing equality 



G{r)d^r 



mm{s,R^f)P{s)ds. 



(37) 



Both sides are proportional to the mean of the length of 
that portion of the contours passing through origin which 
lies within a radius R from the origin [30|. 



Final remark concerns the probability distribution of 
contour length. To this end we investigated the logarith- 
mic diagram of P{s)s^~^ versus s. We have depicted the 
results for the synthetic singular as well as smoothened 
multifractal surfaces for various values of h{q = 2). For 
the smoothened multifractal rough surfaces again the ex- 
ponents follow the behavior of the mono fractal surfaces 
(see Fig. 14). 

In spite of the huge difference between the geonietri- 



VI. CONCLUSION 

In this paper we have studied the contour lines of 
particular multifractal rough surfaces, namely the so- 
called multiplicative hierarchical cascade p model. Uti- 
lizing a stochastic cascade method [J, [SJ] , singular mea- 
sure (original) and smoothened (convolved) multifractal 
rough surfaces with different Hurst exponents were gener- 
ated. The h{q) spectrums of these two dimensional sur- 
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Exponent 


Singular measure 


Smoothened niultifractal 


Xl 


Depends on p— values 


1 
5 


Df 


Identical 


3-H* 
9, 


D{q) 


Depends on q 


Depends on q 


d 


Depends on p— values 


2-H* 


r,-l 


Depends on p— values 


A-2H* 


C 


Depends on p— values 


2-H* 


4j-v-i 


YES 


YES 


Df='-^ 


YES 


YES 



TABLE VI: A summarization of results given in this pa- 
per based on contouring analysis for synthetic singular and 
smoothened multifractal rough surfaces. 



faces were determined by MF-DFA method [34|. Then 
use of two different algorithms we generated the contour 
loops of the systems. Many different geometrical expo- 
nents, such as the fractal dimension of a contour loop, 
D f the fractal dimension of the contour set, d, the cu- 
mulative distributions of perimeters and areas and the 
correlation exponents, xi, were calculated for the singular 
and smoothened multifractal surfaces by use of different 
methods. 

We summarize the most important results given in this 
study as follows. Our results confirmed that, the ex- 
ponent of loop correlation function, x/, for multifractal 
singular measure, depends on p values. On the contrary, 
for multifractal smoothened surfaces, this value behaves 
the same as that of given for mono fractal rough surfaces 
(see Fig.©. 

The scaling exponent of the size of the contours as a func- 
tion of the radius representing fractal dimension, Dy, is 
similar for various singular multifractal rough surfaces. 
But the relation between Df and H* for convolved mul- 
tifractal surfaces is similar to mono fractal surfaces. Nev- 



ertheless the contour loops have multifractal nature (see 
Figs. [TO] and [ni). 

The exponent of cumulative distribution of areas, ^, 
for singular measure has multifractal nature. But for 
smoothened surface, this quantity is controlled by H* 
according to ^ = 2 — H* and is completely independent 
of p values (see Fig. [T3| . 

Consequently, in the case of singular measure surfaces, 
all of the exponents show significant deviations from the 
well-known formulas for the mono- fractal rough surfaces. 
They depend on the generalized Hurst exponents h{q), 
whreas for convolved multifractal surfaces, all geometri- 
cal exponents are controlled by H* according to a mono 
fractal system. We emphasize that interestingly, the hy- 
perscaling relations, namely, -jI- = 77 — 1 and Df — ^^'~^ 



>)-3 



at la confidence interval, are valid for both singular and 
smoothened multifractal rough surfaces(see Table. IIVI 
|V|) . In this system which is labeled by H* , many rel- 
evant properties are controlled by a few relations that 
have been presented for mono-fractal cases. However sin- 
gular and smoothened multifractal surfaces have multi- 
fractal nature but using geometrical analysis, they belong 
to different class which is a non-trivial result. Table IVll 
contains most important results given in this paper. 

Finally, to make more complete this study, it is use- 
ful to extend this approach for various simulated rough 
surfaces by use of different methods and examine their 
hyperscaling relations. In addition, there are some meth- 
ods to distinguish various multiplicative cascade methods 
such as n-point statistics |6C| . 
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